Dynamical Mean-Field Theory for Molecular Electronics: 
Electronic Structure and Transport Properties 
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We present an approach for calculating the electronic structure and transport properties of 
nanoscopic conductors that takes into account the dynamical correlations of strongly interacting 
dr or /-electrons by combining density functional theory calculations with the dynamical mean-field 
theory. While the density functional calculation yields a static mean-field description of the weakly 
interacting electrons, the dynamical mean-field theory explicitly takes into account the dynamical 
correlations of the strongly interacting d- or /-electrons of transition metal atoms. As an example 
we calculate the electronic structure and conductance of Ni nanocontacts between Cu electrodes. 
We find that the dynamical correlations of the Ni 3d-electrons give rise to quasi-particle resonances 
at the Fermi-level in the spectral density. The quasi-particle resonances in turn lead to Fano line- 
shapes in the conductance characteristics of the nanocontacts similar to those measured in recent 
experiments of magnetic nanocontacts. 

PACS numbers: 73.63.-b,73.63.Rt,75.47. Jn 



I. 



INTRODUCTION 



State of the art for calculating the conductance and 
current through atomic- and molecular-size conductors 
consists in combining ab initio electronic structure cal- 
culations on the level of density-functional theory (DFT) 
with the Landauer formalism or non-equilibrium Green's 
function technique^ 2 -. This methodology has been quite 
successful for the theoretical description of e.g. metallic 
nanocontacts 3 - predicting zero-bias conductances that are 
in general in good agreement with experiments even in 
the case of nanocontacts made from magnetic transition 
metals^—. 

Recently, however, nanocontacts made from Fe, Co, 
or Ni have been reported to display Kondo effec t 10 ' 11 . 
This has been infered from the observation of Fano- 
lineshapes 1 ^ in the low-voltage conductance characteris- 
tics similar to those observed in recent STM experiments 
with magnetic adatoms on metal surfaces^—. The ob- 
servation of the Kondo effect in Fe, Co, and Ni nanocon- 
tacts is rather surprising since these materials are strong 
ferromagnets as bulk materials. The Kondo effect, how- 
ever, is usually at odds with ferromagnetism as it results 
from the screening of a local magnetic moment by anti- 
ferromagnetic coupling to the conduction electron s 10 ' 17 . 

DFT based transport calculations of nanoscopic con- 
ductors cannot capture truly many-body effects as the 
Kondo effect that originate from the dynamic correla- 
tions of strongly interacting electrons. Therefore it is 
necessary to extend the existing methodology in order to 
account for dynamic electron correlations in transport ex- 
periments of nanoscopic conductors. Recent approaches 
to include dynamic electron correlations in the ab ini- 
tio description of quantum transport are based on the 
GW approximation (GWA) 1 — or the three-body scatter- 



ing formalism (3BS)~ While the GWA is only suitable 
for weakly correlated systems due to the perturbativc 
treatment of the electron-electron interactions, the 3BS 
is in principle capable of describing more strongly cor- 
related systems as it goes beyond perturbation theory. 
However, the 3BS does not provide a satisfactory solu- 
tion of the Anderson impurity problem since the local 
correlations are not taken into account properly. More re- 
cently, ab initio electronic structure methods on the level 
of DFT or the GW approximatuion have been combined 
with more sophisticated many-body techniques such as 
the One-Crossing Approximation (OCA) or the Numer- 
ical Renormalization Group (NRG) in order to account 
for the Kondo effect in nanoscopic systems containing 
single magnetic atoms^Sr— 

In this paper we develop an approach inspired by the 
success of the Dynamical Mean-Field Theory (DMFT) 
in the treatment of correlated solids 24,25 , to tackle the 
challenges of molecular electronics. DMFT is an ap- 
proach based on the locality of the self-energy and does 
not require Bloch periodicity. It has been applied be- 
fore to strongly spatially inhomogeneous systems such 
as alloys near an Anderson transition 2 ^, surfaces^ 7 - and 
interfaces 2 ^, multilayered heterostructures 29 and cold 
atoms in a trap^ -. Notice, however, that so far all these 
studies were restricted to model Hamiltonians as for ex- 
ample the Hubbard model. 

In order to incorporate realistic aspects of the 
electronic structure we extend the DFT+DMFT 
philosophy 31 ' 32 to the case of nanoscopic conductors. 
Our Molecular DMFT approach takes into account the lo- 
cal dynamical correlations of the strongly interacting 3rf- 
or 4/-shells of the magnetic atoms within a nanoscopic 
conductor such as a molecule or nanocontact which is 
coupled to semi-infinite electrodes while the rest of the 
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system is treated on a static mean-field level by the local 
density approximation (LDA). This approach, can also 
be viewed as an extension of the early work of one of 
us 3 - 3 - which ignored dynamical correlations of open atomic 
shells. In the limiting case where the correlated region 
reduces to a single atom, our approach reduces to our pre- 
vious work which treated a single magnetic impurity in a 
metallic nanocontact^. When the device region contains 
several correlated atoms in close proximity, the Molecular 
DMFT treatment is essential since the effective bath of 
each correlated atoms obeys a self-consistency condition 
which involves the whole device region. 

Our work is similar in spirit to the recently presented 
dynamical vertex approximation for nanoscopic systems 
(nano-DrA) by A. Valli et al£^ This approach can in 
principle treat non-local self-energy effects. In a re- 
lated work, Florens introduced a nano-DMFT approach 
whereby a correlated system is approximated by embed- 
ding it in a model geometry having a tree like structured 
However, both method a 34 i 35 have been implemented in 
the context of the Hubbard model while we demonstrate 
that the existing Molecular DMFT technology can treat 
a realistic nanocontact. 

This paper is organized as follows. In Sec. |H] wc de- 
scribe the Molecular DMFT method for nanoscopic con- 
ductors. In Sec. IIIII we apply the method to small Ni 
nanocontacts connected to Cu nanowires, and discuss the 
results. Finally, we conclude this paper with a general 
discussion of the method and of the results in Sec. IIVI 



II. METHOD 

We consider a nanoscopic conductor bridging two semi- 
infinite metal wires as shown schematically in Fig. [TJb). 
As indicated, the nanoscopic conductor contains mag- 
netic atoms that give rise to dynamic electron correla- 
tions due to the strongly interacting 3d-electrons. The 
nanoscopic conductor could be for example a molecule, a 
nanowire, a nano-cluster, or simply an atomic-size con- 
striction in the metal wire. 

In order to describe the dynamic correlations that arise 
from the 3d-shells of the magnetic atoms we adapt the 
LDA+DMFT method to the case of nanoscopic conduc- 
tors. To this end it is convenient to divide the system 
into three parts as shown in the upper right panel of Fig. 
[TJ Two semi-infinite metallic leads L and R, and the 
central device region (D) which contains the nanoscopic 
conductor and the magnetic atoms with the strongly in- 
teracting 3<i-shells, as well as a sufficient part of the two 
metallic leads so that the electronic structure of the leads 
has relaxed to that of bulk (i.e. infinite) electrodes. The 
correlated subspace (C) that will be described on the 
level of DMFT is thus given by the direct sum of the 3d- 
subspaces of all magnetic atoms (M) in the device region. 

The effective one-body Hamiltonians of the device re- 
gion and leads are obtained from DFT calculations on 
the level of LDA. Here we use the supercell approach 
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FIG. 1: (Color online) (a) Molecular DMFT self-consistency 
cycle applied to nanoscopic conductors, (b) Division of the 
system into left (L) and right (R) electrodes, and the central 
device region (D) that hosts the strongly correlated subspace 
(C) consisting of the strongly correlated d-orbitals of the mag- 
netic atoms shown as red (dark grey) circles. 



(see App. [A| to obtain the effective Kohn-Sham (KS) 
Hamiltonians of each part of the system prior to the dy- 
namical treatment of the impurity d-shell and the trans- 
port calculations. The electronic structure of the device 
region is calculated with the Crystal06 ab initio elec- 
tronic structure program for periodic systems 3 -^ by defin- 
ing a one-dimensional periodic system consisting of the 
device region as the unit cell. 

The device Hamiltonian Hd is then obtained from the 
converged KS Hamiltonian of the unit cell of the peri- 
odic system. In the same way, the unit cell Hamiltoni- 
ans H"^ R and hoppings Vl/r between unit cells of the 
left and right leads can be extracted from calculations of 
infinite nanowires with finite width since the electronic 
structure in the semi-infinite leads has relaxed to that of 
an infinite nanowire. 

The strong electron correlations of the 3d-shells of the 
magnetic atoms are captured by adding a Hubbard-type 
interaction term to the one-body Hamiltonian within the 
correlated subspace d of each of the magnetic atoms i: 
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where d iacr (di a<J ) creates (annihilates) an electron with 



spin a in the 3<i-orbital a on atom i. c 

r aiOt2PlP2 

matrix elements of the effective (i.e. screened) Coulomb 
interaction of the 3d-electrons which is smaller than the 
bare Coulomb interaction due to the screening by the 
conduction electrons. Here we take a model interaction 
taking into account only the direct Coulomb repulsion 
U between electrons (i.e. U a p a @ = U) and the Hund's 
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rule coupling J (i.e. U a fjp a = J for a ^ (3). For 3d 
transition-metal elements in bulk materials the repulsion 
U is typically around 2-3 eV and J is around 0.9 eV.— 
However, due to the lower coordination of the atoms in 
the contact region or molecule the screening of the di- 
rect interaction should be somewhat reduced compared 
to bulk. Here we take U = 5 eV and J = 0.9 eV as in 
our previous work^S 

The Coulomb interaction within the correlated 3d sub- 
spaces of the magnetic atoms has already been taken 
into account on a static mean-field level in the effec- 
tive KS Hamiltonian of the device. Therefore the KS 
Hamiltonian within each 3d subspace has to be cor- 
rected by a double-counting correction term, i.e. = 

Pj^H^P^ - where is the projection onto 
the 3d-subspace of atom i.— Here we use the standard 
expression 

H« = [U(N® -I)- \j{Nf 1)] x pW (2) 

where N^' is the occupation of the 3(i-shell of atom i. 

The central quantity is the Green's function of the de- 
vice region: 

G D = ((w + M )S D - H D + H dc - S c - S L - SrT 1 

(3) 

where Sc is the electronic self-energy that describes the 
dynamic electron correlations of the electrons within the 
correlated subspace C, and Hd c is the double-counting 
correction within the entire correlated subspace C, i.e., 
Hdc = YlieM Hdc- So is the overlap matrix taking into 
account the non-orthogonality of basis set4S ^ is the 
chemical potential, Sl and Sr. are the so-called lead 
self-energies^ which describe the coupling of the device 
to the semi-infinite leads L and R, respectively. These 
can be calculated from the effective one-body Hamiltoni- 
ans of the leads by iteratively solving the Dyson equation 
for the lead self-energies, eq. (|Af [) . 

The central assumption of DMFT is that the inter-site 
correlations, i.e. the correlations between electrons lo- 
cated on different atoms can be neglected. In that case 
the electron correlation self-energy Sc becomes block- 
diagonal, and each block corresponds to the self-energy 
Sy of the correlated 3d subspace of a magnetic atom 

i: = P^ScP^. The self-energies and hence 
the overall self-energy Sc can now be determined by 
mapping onto a generalized Anderson impurity problem 
for each correlated subspace, described by the following 
Green's function: 

= (" + M - H« - EW( W ) _ A (0 (w) ) 1 (4) 

where we have introduced the so-called hybridization 
function Aj' which describes the coupling of the cor- 
related subspace with the rest of the system. The hy- 
bridization function is determined by the DMFT self- 



consistency condition (DMFT-SCC): 

g«(u;)=P«G D (u;)P« (5) 
It follows that the hybridization function is given by: 

A«H = w + M - H« - [gf (a;)]- 1 - (6) 

The equations (0O define the Molecular DMFT self- 
consistency cycle for the calculation of the self-energies 
: One starts with the effective one-body Hamiltonian 
Hd obtained from the LDA calculation and an initial 
guess for the local correlation self-energies Ejj (usually 
zero). This allows one to calculate the device Green's 
function, eq. ([3]), and consequently the projection gl . 
Thus one obtains the Hybridization functions Aj (cj) 

(i) 

which together with the on-site energy levels and 
the Coulomb interaction U, J defines the Anderson impu- 
rity problem which can be solved by an impurity solver, 
and one obtains a new self-energy Yl d . The DMFT self- 
consistency cycle is illustrated in Fig. Qla). 

Solving the generalized Anderson impurity problem is 
a difficult task, and at present there is no universal im- 
purity solver that works efficiently and accurately in all 
parameter regimes. Here we make use of impurity solvers 
based on an expansion in the hybridization strength given 

by (w) around the atomic limit. The starting point is 
an exact diagonalization of the (isolated) impurity sub- 
space i.e., the 3d-shell of the magnetic atom given by the 

interacting Hamiltonian H d +Hjj ■ The hybridization of 
the impurity subspace with the rest of the system (given 
by the hybridization function Aj (w)) is then treated 
perturbatively. 

The so-called Non-Crossing Approximation 42 (NCA) 
is a self-consistent perturbation expansion to lowest or- 
der in the hybridization strength. NCA only takes into 
account bubble diagrams describing hopping processes 
where an electron or hole hops into the impurity at some 
time and then out at a later time (see Fig. [5] in App. 
HJ. The OCA 4 ^ improves on the NCA by taking into 
account second order diagrams where two additional elec- 
trons (holes) are accommodated on the impurity at the 
same time as shown in Fig. [9] OCA improves consider- 
ably many of the shortcomings of NCA 42 : It substantially 
improves the width of the Kondo peak and hence the 
Kondo temperature which now are only slightly under- 
estimated. It also corrects the asymmetry of the Kondo 
peak. For very low temperatures (T <C Tk), however, 
the height of the Kondo peak is still overestimated, and 
the Fermi liquid behavior at zero temperature is not 
recovered^ 

Hence, OCA is a reasonable approximation for solving 
the generalized impurity problem as long as the tempera- 
tures are not too low (i.e. more than one order of magni- 
tude below Tk)- In App. |B]we give a brief introduction 
to the NCA and OCA impurity solvers. A detailed de- 
scription of the NCA and OCA methods can be found 
e.g. in Refs. [I3|3l]|I3fil 
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FIG. 2: (Color online) (a) Ni dimer between the tips of two 
Cu nanowires grown in the 001 direction of bulk Cu. (b) 
Ni nanocontact consisting of 10 atoms between two Cu 001 
nanowires. Ni atoms are shown in light grey (light grey) and 
Cu atoms are shown in red (dark grey). 
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The current through a strongly interacting region can 
be calculated exactly by the Meir-Wingreen formula 46 . 
But in order to apply the Meir-Wingreen result one 
has to solve the impurity problem out of equilibrium 
which is a difficult task that has only been accom- 
plished very recently and only in the context of model 
Hamiltonian a 47 i 48 . However, Meir and Wingreen also 
showed that for low temperatures and small bias volt- 
ages the Meir-Wingreen expression is well approximated 
by the much simpler Landauer formula^ 

W = f x £ V dcjTiu) (7) 

where T(u>) is the Landauer transmission function and 
where we have assumed an asymmetric voltage drop V 
about the device region. Thus the conductance is simply 
given by the Landauer transmission function: 

31 2e 2 
G(V) = ^(V) = — xT(eV) (8) 

The latter can be calculated from the (equilibrium) de- 
vice Green's function: 

T(w) = Tr[r L HGfy w )r R HG D H] (9) 

where Tl/r are the so-called coupling matrices which 
describe the coupling to the leads, and can be calculated 
from the lead self-energies by r L /R, = i(S l /r, — , R ). 



III. RESULTS AND DISCUSSION 

In the following we demonstrate the above devel- 
oped Molecular DMFT method for the two idealized Ni 
nanocontacts with Cu nanowires as electrodes shown in 
Fig. [5] Ni nanocontacts have been a subject of intense 
research in the last decade as prospective ingredients for 
nanoscale spintronics devices^ - — i 50 ' 51 Here we consider 
the paramagnetic phase, i.e. the self-energies and the 
Green's functions are spin-degenerate. Breaking of the 
spin-symmetry by an external magnetic field or by spin- 
polarized electrodes is not taken into account. 
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FIG. 3: Imaginary part of hybridization functions for differ- 
ent Ni 3d-orbitals in the dimer geometry shown in Fig. [2{a) in 
the first step (a) and in the last step (b) of the self-consistent 
Molecular DMFT calculation at low temperature (T = 12 K). 



A. Ni dimer between Cu nanowires 

First, we consider a dimer of Ni atoms between two 
semi-infinite Cu nanowire electrodes as shown in Fig. 
The nanowires are grown in the 001 direction of bulk Cu. 
The distance between the two Ni tip atoms is 2.4 A, all 
other distances are those of bulk Cu. For the sake of sim- 
plicity, we have not relaxed the atomic positions. Due to 
the highly idealized geometry both Ni atoms are equiv- 
alent. Hence in each step of the Molecular DMFT self- 
consistency we only have to solve the impurity problem 
once. 

In spite of the highly symmetric situation the Ni 3d- 
orbitals split into four different symmetry groups. As 
can be seen from Fig. EHa) which shows the hybridiza- 
tion function in the first step of the self-consistent DMFT 
procedure (where the self-energy Sd(w) is zero), the hy- 
bridization functions are quite different for each of the 
four symmetry groups. This situation is different from 
that of the corresponding bulk systems where the hy- 
bridization functions for each correlated orbital are usu- 
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FIG. 4: Ni 3d-spectral function for the Ni dimer between 
Cu nanowires shown in Fig. [2ja) calculated with Molecular 
DMFT at two different temperatures on the one hand and 
calculated with LSDA on the other hand. 



ally very similar due to the highly isotropic closed packed 
crystal structures. In the geometry considered here, 
the doubly-degenerate 2>d xz - and id yz - orbitals have the 
strongest hybridization around the Fermi-level with the 
rest of the system. Also the 3d 3z 2_ r 2 -orbital has an ap- 
preciable hybridization. The hybridization of the 3d xy - 
and 3d x 2_ y 2 -orbitals on the other hand are smaller by at 
least one order of magnitude. 

Figure^b) shows the hybridization function after self- 
consistency has been reached in the DMFT calculation. 
Now the converged self-energy T,d(u) is non-zero in gen- 
eral. We can see that the DMFT self-consistency has 
a considerable effect on the hybridization function for 
most of the Ni 3<i-orbitals. For example, the hybridiza- 
tion of the degenerate 3d xz - and Sii^-orbitals around the 
Fermi level is strongly increased. Moreover a sharp peak 
appears right at the Fermi level. Additional features 
also arise in the hybridization functions of the other 3d- 
orbitals. The additional features in the 3d-hybridization 
functions of one Ni atom stem from the corresponding 
features (induced by the on-site interactions) in the 3d- 
spectral density of the other Ni atom. 

In Fig. @]we show the temperature dependence of the 
Ni 3d-spectrum calculated with Molecular DMFT and 
compare them to the spectrum calculated with LSDA. 
The DMFT spectra are somewhat smother than the 
LSDA spectra due to the finite lifetime broadening of 
the single-particle states by the electron-electron interac- 
tions. Most importantly, the Molecular DMFT spectrum 
shows a strong temperature dependence which cannot be 
captured by a static mean-field treatment like LSDA. For 
example, at low temperatures a small peak forms right at 
the Fermi level in the Molecular DMFT spectrum. This 
is a quasi-particle peak that originates from the two de- 
generate ?>d xz - and 3d yz -orbitals as can be seen from Fig. 
[SJa). Strictly speaking, it is not a Kondo peak since the 
occupation of the two orbitals is around 3.65, and hence 




FIG. 5: (a) Spectral function of Ni 3d xz - and Sd^z-orbitals 
near the Fermi level for different temperatures, (b) Correlated 
transmission function calculated with Molecular DMFT at 
different temperatures compared to LSDA transmission func- 
tion near the Fermi level. 



these orbitals are in the so-called mixed-valence regime 
where in addition to the spin-fluctuations also charge 
fluctuations take place (see e.g. the book by Hewson 
[l7l for a detailed discussion). However, we would like to 
emphasize here that the corresponding magnetic moment 
is nevertheless screened by the conduction electrons. 

We can find the weight Z and the linewidth V of the 
quasi-particle by fitting the peak to a Lorentzian. Fig. 
[5{a) suggests that the 3(4,2-, 3dj, z -spectral function near 
the Fermi level can be fitted to a weighted Lorentzian 
plus a linear function: 



A d (u) 



T/2 



7TCJ 2 +r 2 /4 



(10) 



We find a very tiny quasi-particle weight Z of less than 
0.1%, and a width T corresponding to a temperature of 
about 130 K. This is the critical temperature below which 
the quasi-particle can be observed (and which in the case 
of Kondo effect is called Kondo temperature). Note that 
the width of the quasi-particle is somewhat enhanced by 
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the DMFT self-consistency compared to the case without 
DMFT (not shown). 

The quasi-particle peak leads to a corresponding Fano 
feature in the transmission function as can be seen in Fig. 
[5jb) which shows the transmission function for small en- 
ergies around the Fermi level. Hence our calculations 
show that the low bias conductance of the Ni dimer be- 
tween Cu electrodes features a Fano-lineshape due to a 
quasi-particle peak in the Ni 3d-spectral function at low 
temperatures. 

Table Q] lists the individual occupations of all Ni 3d- 
orbitals and the corresponding effective energy levels that 
are obtained from the KS energy levels plus the real part 
of the hybridization function at zero frequency. We can 
see that all orbitals apart from the d^-orbital have mixed 
valences and are almost full. The d xy -orbital on the other 
hand is basically half-filled, hence this orbital would be 
a candidate for a true Kondo effect, i.e. screening of the 
magnetic moment by spin-fluctuations only. However, 
the hybridization of this orbital is very low, as can be 
seen from Fig. [3] Hence the Kondo temperature for this 
orbital is very low so that the orbital does not enter the 
Kondo regime at the temperatures considered here but 
is in the local moment regime. 



B. Ni nanocontact between Cu electrodes 

Now we turn to the slightly more complicated case 
of the Ni nanocontact consisting of 10 atoms between 
two Cu nanowires as shown in Fig. [2jb). As before the 
nanowires are grown in the 001 direction of bulk Cu. 
The distance between the two Ni tip atoms is 2.4 A 
while all other distances are those of bulk Cu. As be- 
fore the two Ni tip atoms are both equivalent. On the 
other hand, the eight outer atoms of the Ni nanocontact 
are not equivalent with the tip atoms, but are equivalent 
among themselves. Hence we have to solve two different 
impurity problems in each step of the Molecular DMFT 
self-consistency cycle. 

In Fig. |S1 we show the 3d-spectral function of (a) the 
base atoms and (b) the tip atoms of the Ni nanocontact 
calculated on the one hand with Molecular DMFT at dif- 
ferent temperatures, and on the other hand with DFT on 
the level of the LDA. We see that in both cases the Molec- 
ular DMFT spectral densities are quite different from the 
LSDA ones. Both LSDA spectra feature a strong peak 
above the Fermi level for the minority electrons that is 
absent in the Molecular DMFT spectra. Moreover, near 
the Fermi level, the Molecular DMFT spectrum of the tip 
atoms is strongly temperature dependent due to the for- 
mation of a quasi-particle peak at the Fermi level in the 
degenerate 3d xz - and 3d yz -levels for low temperatures. 
Fig. 0(a) shows a close-up of the formation of the quasi- 
particle peak in the spectral density of the 3d xz - and 
3cij, z -levels. Fitting to a Lorentzian plus linear function, 
cq. (|10[) . we now find a slightly increased quasi-particle 
weight Z of 0.2% and also a somewhat increased width T 




FIG. 6: (a) 3d-spectral function of a Ni atom at the base 
of one of the pyramids of the Ni nanocontact between Cu 
nanowires shown in Fig. EJb) calculated on the one hand 
with Molecular DMFT at two different temperatures and with 
LSDA on the other hand, (b) Same as (a) but for the tip 
atoms of the Ni nanocontact. 



corresponding to a critical temperature of 220 K as com- 
pared to the case of the Ni dimer. The increased width T 
and quasi-particle weight Z is due to the increased imag- 
inary part of the hybridization function (not shown) of 
the tip atoms near the Fermi level which stems from the 
3d-spectral density of the base atoms of the Ni pyramid. 
In the dimer case, the base atoms of the pyramids are Cu 
atoms where the 3d-spectral density near the Fermi level 
is negligible. 

The occupation of the degenerate 3d xz - and 3d yz - 
orbitals is 3.73. Hence these orbitals are in the mixed- 
valence regime rather than the Kondo regime as in the 
case of the Ni dimer. Therefore the quasi-particle peak is 
strictly speaking not a Kondo peak. However, as was said 
before, the magnetic moment of these orbitals is never- 
theless screened in the mixed valence regime by the spin- 
and charge-fluctuations. 

Figure [TJb) shows the transmission function near the 
Fermi level calculated on the one hand with the Molecu- 
lar DMFT method at different temperatures and on the 
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FIG. 7: (a) Spectral function of 3d xz - and 3d yz - orbitals of 
Ni tip atoms near the Fermi level for different temperatures, 
(b) Correlated transmission function at different temperatures 
compared to LSDA transmission function for energies around 
the Fermi level. 

other and with DFT on the level of the LDA. As in the 
case of the dimer, the formation of the quasi-particle peak 
in the 3d xz - and 3<ij, 2 -orbitals at low temperatures leads 
to a temperature dependent Fano lineshape in the Molec- 
ular DMFT transmission function. The LSDA transmis- 
sion function on the other hand does not show such a 
behavior. On the contrary, it is rather flat and feature- 
less at this energy scale. Also note that the LSDA trans- 
mission is considerably higher than the correlated trans- 
mission calculated with Molecular DMFT. This is due to 
the correlations shifting a considerable part of the spec- 
tral weight of the 3c?-orbitals away from the Fermi level. 
Consequently, the overall contribution of the 3<i-orbitals 
to the transmission is higher in the case of the LSDA 
calculation than in the Molecular DMFT calculation. 

The transmission as calculated with the Molecular 
DMFT approach, although lower than the LSDA one, is 
still compatible with the broad peak roughly between 1 
and 1.7xGo in the experimentally measured conductance 
histograms of Ni nanocontacts£ 

In Tab. [I] the orbital occupations and effective energy- 
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-4.8121 


1.0096 


-5.1749 


xy 


0.9943 


-4.7410 


1.9578 


-5.2706 



TABLE I: Orbital occupations rid (at low temperature) and 
effective energy levels ed = ed + ReAd(0) of the 3d-orbitals of 
the Ni tip atoms in the case of the Ni dimer geometry and 
the Ni nanocontact consisting of 10 Ni atoms, ed denotes the 
KS energy levels before double-counting correction. 

levels of individual Ni 3d-orbitals of a tip atom of the 
Ni nanocontact are shown. The most striking difference 
with the Ni dimer (shown in the same table) is that now 
the d x 2_ y 2 -orbital is the highest energy orbital with half- 
filling instead of the d xy - orbital. This is a consequence of 
the different environment of the tip atom in the 10-atom 
nanocontact geometry compared to the dimer geometry. 
The presence of the 3d-orbitals of the Ni atoms at the 
base of each pyramid near the Fermi levels change the hy- 
bridization functions of the tip atoms accordingly. Note 
that here the DMFT self-consistency is essential since 
without the self-consistency the hybridization functions 
would be very similar to the dimer case. 

IV. SUMMARY AND CONCLUSIONS 

We have developed a method for calculating the elec- 
tronic structure and transport properties of nanoscopic 
conductors that explicitly takes into account local dy- 
namical correlations originating from strong electron- 
electron interactions. Our method extends the estab- 
lished DFT based ab initio transport methodology for 
nanoscopic conductors to include dynamic electron cor- 
relations originating from the strongly interacting 3d- 
electrons of the transition metal atoms. This is achieved 
by combining the DFT electronic structure calculations 
of the nanoscopic conductor with a DMFT description 
of the strongly interacting 3<i-electrons in the device re- 
gion. We thus obtain the correlated Green's function of 
the nanoscopic conductor which allows to calculate the 
electronic structure and the corresponding conductance 
in the low-bias voltage regime. 

We have demonstrated the method for two model 
systems, namely Ni nanocontacts consisting of several 
atoms and connected to Cu leads. We find that the 
dynamic correlations of the strongly interacting Ni 3d- 
electrons give rise to strongly temperature-dependent 
spectra due to the formation of a quasi-particle peak 
at low temperatures. The quasi-particle peak gives 
rise to a temperature-dependent Fano-type lineshape 
in the low-bias conductance characteristics similar to 
those measured in recent experiments with ferromagnetic 
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nanocontactsAi Moreover, the critical temperatures of 
120 K and 220 K for the formation of the quasi-particle 
peak is in quite good agreement with the broad distri- 
bution of Kondo temperatures around the average tem- 
perature of 250 K extracted from the Fano lineshapes in 
the low-bias conductance measurements of Ni nanocon- 
tacts in that same experiment. Note that a quasi-particle 
peak at the Fermi level can in principle also be obtained 
in the GW approximation^ However, in order to cap- 
ture the quasi-atomic features characteristic of the strong 
correlation regime such as Hubbard bands or satellites to- 
gether with the concomitant renormalization of the quasi- 
particle, a non-perturbative treamtent of the local part 
of the Coulomb interaction such as the Molecular DMFT 
method presented here is necessary. 

The quasi-particle peak obtained here, is strictly 
speaking, not a Kondo peak since the system is in the 
so-called mixed- valence regime where charge fluctuations 
take place in addition to the spin-fluctuations that lead 
to the Kondo effect. This hints at the possibility that the 
origin of the Fano lineshapes in the low-bias conductance 
of Ni nanocontacts measured experimentally need not al- 
ways be the Kondo effect. However, we would like to 
stress that also in the mixed- valence regime the magnetic 
moment of the corresponding orbitals would be screened. 
But to draw further conclusions in that matter, more re- 
alistic calculations are necessary taking into account the 
ferromagnetic leads, and sampling over different contact 
geometries. 

We have illustrated the Molecular DMFT method for 
the case of simple nanocontacts containing several tran- 
sition metal atoms, but the Molecular DMFT approach 
is very general and can be applied to many systems of 
great theoretical and practical interest. For example, it 
can be used to treat large molecules in which one can iso- 
late small clusters of correlated elements as for example 
in the fuel cell materials of Tard et. alM- 

The Molecular DMFT method allows to explicitly in- 
corporate strong dynamical correlations within the estab- 
lished DFT based transport methodology for nanoscopic 
conductors. Our calculations show that dynamical corre- 
lations originating from the strongly interacting shells of 
magnetic atoms can alter the electronic structure and 
transport properties of nanoscopic conductors signifi- 
cantly. 
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FIG. 8: Illustration of the supercell approach to calculate 
the electronic structure of the device and of the leads: (a) 
One-dimensional periodic system to calculate the electronic 
structure of the device region. (b,c) Infinite nanowires to 
calculate the electronic structure of the left (L) and right (R) 
semi-infinite leads, (d) Sketch of the setup of the physical 
system: The device region (D) is suspended between two semi- 
infinite leads L and R. 



Appendix A: Details of the supercell approach 

In order to generate the effective one-body Hamilto- 
nians of the device and leads the supercell approach is 
used. The electronic structure of the device region is cal- 
culated with the Crystal06 ab initial electronic struc- 
ture program for periodic systems by definining a one- 
dimensional periodic system consisting of the device re- 
gion as the unit cell, as shown in Fig. Ufa). It is cru- 
cial that the device part D contains a sufficient part of 
the nanowire electrodes so that the two leads L and R 
are far enough away from the scattering region, and the 
electronic structure in the leads has relaxed to that of a 
bulk (i.e. infinite) nanowire. In that case the electronic 
structure of the periodic system build from a periodic 
repetition of the device region is the same as the elec- 
tronic structure of the device between two semi-infinite 
nanowires as can be seen from Fig. [8] Thus the device 
Hamiltonian Hd can be obtained from the converged KS 
Hamiltonian of the unit cell of the periodic system. 

In the same way the unit cell Hamiltonians H^ R and 
hoppings Vl/r between unit cells of the left and right 
leads are extracted from periodic calculations of infinite 
nanowires of finite width (see Fig. EtR,c)). Again it is 
crucial that the device region contains enough bulk elec- 
trode material so that the electronic structure in the elec- 
trodes is that of bulk nanowires. The lead self-energies 
Sl, which describe the coupling of the device region 
D to the semi-infinite nanowires L and R in the situa- 
tion depicted in Fig. |8{d) can now be calculated by the 
following Dyson equation: 

S l /r(w) = (V l /r - wS L/R ) x (Al) 
x(u,S L/R - H L/R - Sl/rM)- 1 (V£ /R - wS L/R ) 

where S L , R and Sl/r are the overlap matrices taking into 
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FIG. 9: Diagrams for pseudo-particle self-energies in the 
NCA (first row) and OCA (second and following rows) for 
some pseudo-particle m. 

account the non-orthogonality of the basis set within the 
unit cell and between unit cells, respectively. 

By this procedure we have connected the device re- 
gion D with two semi-infinite nanowires that have the 
electronic structure of bulk i.e. infinite nanowires. The 
supercell approach and the so-called partitioning tech- 
nique used here to obtain the Green's function of a part 
of a system are discussed in more detail in the literature 
(See e.g. refs. HHHI). 

Appendix B: The NCA and OCA impurity solvers 

The general multi-orbital Anderson impurity model 
can be written in the following form: 

H = d\d a + - ^ U a ^gd\d^d^ds (Bl) 

hva. kv 

where in order to keep the notation simple we have com- 
bined the spin- and orbital degrees of freedom into one 
index for each impurity level a and each band v. 

The NCA and the OCA both solve the Anderson im- 
purity model by expansion in the hybridization strength 
around the atomic limit. The starting point is an exact 
diagonalization of the impurity subspace (for example the 
3tf-shell of a magnetic atom) including the Hubbard-type 
interaction term: 

Hd = X! £q ^ a da + ^ X! U a fasd\dpd 7 ds 

d ^ H E ™ H (B2) 

m 



where \m) are the many-body eigenstates of T-Ld and E m 
the respective eigen energies. 

One now introduces auxiliary fields a m , (called 
pseudo-particles) such that each impurity state is rep- 
resented by a corresponding pseudo-particle: 

al |PPV) = \m) (B3) 

where |PPV) is the pseudo-particle vacuum. The com- 
pleteness of the impurity eigenstates imposes the follow- 
ing constraint: 

Q = Y j a[ n a m = \ (B4) 

m 

The physical electron operators ct^ can now be expressed 
by the pseudo-particle (PP) operators: 

d\ = ^(F a t) nm at a m (B5) 

n,m 

where (F a ') nm = (nj rfj,, \m) are the matrix elements of 
the impurity-electron creation operator. For later con- 
venience we also define the corresponding matrix ele- 
ments of the impurity-electron destruction operator as: 
(F a ) nm = (n\d a \m). The anti-commutation rules for 
the physical electron operators then require that the PP 
cijii is a boson (fermion) if the corresponding state |m) 
contains an even (odd) number of electrons. 

In the PP representation we can now rewrite the 
Hamiltonian of the generalized Anderson impurity model 
as follows: 

H = EmO^am + e k yc\ v c kv + X(Q - 1) 

m kv 

+ Yl ( v ku,A^li Fa )nma n + H.c)j (B6) 

where we have included the constraint Q = 1 into the 
Hamiltonian. The corresponding Lagrange multiplier A 
can be interpreted as a (negative) chemical potential for 
the PPs. 

In the PP picture, the hybridization with the bath elec- 
trons given by the last term in eq. (|B6[) becomes now the 
interaction for the PPs. Because of the fermionic and 
bosonic commutation rules for the PPs, one can now de- 
velop a diagrammatic perturbation expansion in the PP 
interaction. The PP propagators can be written as 

G m {L>) = (uj-\-E m -X m (oj))- 1 (B7) 

where S m is the PP self-energy describing the dynamic 
interaction of the PP m with the other PPs. 

The NCA consists in taking into account the diagrams 
shown in the first row of Fig. [5] for some PP m. The 
NCA diagrams describe processes where a single elec- 
tron (hole) jumps from the bath to the impurity and 
back thereby temporarily creating a PP with N+l (N- 
1) electrons. The NCA equations correspond to a self- 
consistent perturbation expansion to lowest order in the 



10 



hybridization function A a (uj) = J2k,v^kv,a^ku,a- Since 
the fermionic self-energies depend on the dressed bosonic 
propagators, and vice versa, the NCA equations have to 
be solved self-consistently. Once the NCA equations are 
solved the physical quantities can be calculated from the 
PP self-energies. 

The OCA takes into account second order diagrams 



where two bath electron lines cross as shown in the 
last four rows of Fig. [§J The self-energies for the PPs 
again depend on the full propagators of other PPs, and 
hence the OCA equations also have to be solved self- 
consistently. Further details of the NCA and OCA im- 
purity solver can be found e.g. in Refs. Il7ll3lll42| - |45l 
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